Theory of field induced spin reorientation transition in thin Heisenberg films 
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We consider the spin reorientation transition in a ferromagnetic Heisenberg monolayer with a 
second order single ion anisotropy as a function of temperature and external field. Up to now 
analytical methods give satisfying results only for the special case that the external field is aligned 
parallel to the easy axis of the crystal. We propose a theory based on a generalization of the Callen 
decoupling, which can be used for arbritrary direction of the external field. Excellent agreement 
between our results and Quantum Monte Carlo data is found for the field induced reorientation at 
finite temperatures. Additionally, we discuss the temperature dependence of the transition in detail. 



I. INTRODUCTION 

Since the discovery of the Giant Magneto Resistance 
(GMR) effect 1989i there has been enormous interest 
and research activity in the field of thin magnetic films. 
The magnetic anisotropy, merely a small perturbation 
in a bulk ferromagnet, gets strikingly important in thin 
film systems. Here the anisotropy is not only a neces- 
sary precondition of spontaneous ferromagnetism^, but 
it determines many system properties as, e.g., the de- 
pendence of the magnetization vector or of the spin wave 
excitation spectrum on an applied magnetic field. Ad- 
ditionally, the anisotropy energy is of the same order of 
magnitude as the Interlayer Exchange Coupling (IEC) 3 , 
which is intimately connected with the GMR effect. Thus 
an investigation of these effects has to take the magnetic 
anisotropy carefully into account. 

There is an important phenomenon in thin ferromag- 
netic films which is closely connected to the magnetic 
anisotropy: the magnetic reorientation transition. This 
term denotes a rotation of the magnetization from the 
film normal into the plane or vice versa as a function of 
temperature, film thickness, or external magnetic field. 
The transition can be understood as a result of compet- 
ing forces that favor different directions of the magne- 
tization as e.g. spin orbit coupling, dipolar interaction, 
and an external magnetic field 4 . It can be described us- 
ing a Heisenberg model in film geometry, with the usual 
Heisenberg exchange interaction, an external field, and 
one or more anisotropy terms. 
For the simplest of these models, 
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JijSiSj — BgSi 



(1) 



consisting only of the exchange term and an external 
field, there are very accurate approximation schemes 
available. It was shown e.g. in Rcf. 5 by comparison with 
QMC calculations, that the RPA decoupling 7 yields even 
quantitative results for the magnetization as a function 
of temperature. 

Turning to anisotropy contributions, the spin-orbit cou- 
pling induced anisotropy is usually modelled by a single 
ion anisotropy 



which is of second order for systems with tetragonal sym- 
metry. For film systems the z-axis is perpendicular to 
the film plane. If K 2 is positive the easy axis of the mag- 
netization is the z-axis, for negative Kq, this is a hard 
direction. The RPA fails badly if applied to a local term 
as described in Eq. Thus RPA can not be used to 

solve the whole model 



H = H 1 +H 2 . 



(3) 



Ho = -Kn V S iz S, 



(2) 



In Rcf. 8 an approximation for this model is proposed, 
which is based on a combination of the RPA approxima- 
tion for the nonlocal terms QJ and an Anderson-Callen 
(A.C.) decoupling^ for the local anisotropy contribution 
@. This theory gives good results 6 for the magneti- 
zation if the anisotropy constant K 2 is much smaller 
than the exchange coupling J (K 2 < 0.01J) and if the 
external field is applied parallel to the z-axis while K 2 
has to be positive. The first condition is not a serious 
restriction, since in reality the anisotropy constants 
are indeed much smaller than the Heisenberg exchange 
interaction. Furthermore this restriction can be relieved 
by an alternative theory 11 . The important restriction 
is given by the second condition. The described limit 
is a very special one, where both, the anisotropy as 
well as the external field favor a alignment of the 
magnetization parallel to the z-axis. Thus there are 
no " competing forces" , and no magnetic reorientation 
transition occurs in this limit, which we want to refer 
to as "parallel limit" in the following. However, if the 
external field is not applied parallel to the z-axis and 
the magnetization is consequently rotated out of z, the 
approximation described in Ref. 8 looses its accuracy 
and becomes unacceptable for a quantitative description 
of the reorientation transition. This was shown in Rcf. 6 
by comparison with QMC calculations. Actually, to our 
knowledge there is no reliable model theory available, 
which can treat the model for arbritrary directions 
of the external field or a negative anisotropy constant 
K 2 . 

However, such a model theory is highly desirable. It 
can be used to investigate quantitatively the magnetic 
reorientation transition in all systems dominated by 
second order lattice anisotropies. Numerical methods, as 
QMC calculations, are only applicable for the monolayer, 
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but a finite number of layers is crucial to study the 
interplay between surface and bulk anisotropies (see e.g. 
Ref. 4). 

A priori, it is not clear, that such a model theory exists. 
The full model @ is much more complex then the 
"parallel limit" case. The reason is that in that special 
case the total spin is a conserved quantity i.e. the total 
magnetization Sf commutes with the Hamiltonian. 
This property simplifies the calculations considerably 
but it is not present in the general model ©. 
In this paper we want to show that nevertheless a well 
funded model theory can be formulated and that it is as 
accurate in the general case as in the "parallel limit". 
Comparing our results with the QMC data of reference 6 
we will find excellent quantitative agreement for the 
magnetization and its components, which allows for a 
high quality description of the magnetic reorientation 
transition. The theory proposed in Ref. 8 is recovered 
for the special case of the " parallel limit" . 
In this paper we want to introduce the new theory 
and evaluate it by comparison with available QMC 
results. Since the latter are results for the monolayer, 
we will exclusively treat the case of a monolayer during 
this paper. A generalization to a multilayer system is 
straightforward. 



II. THEORY 

We will assume tetragonal symmetry in the following. 
The xy-plane is the film plane and thus the x— and the 
y— direction are equivalent. That is why we can confine 
the external field and the magnetization to the zcc-plane 
without loss of generality. We will assume nearest neigh- 
bor coupling (Jij — J for nearest neighbors and J = 
elsewhere) and for the explicit calculations a quadratic 
lattice. Let us first outline the main points of our theory. 
The aim is to calculate the angle and the norm of the 
total magnetization of the model Q. 

1. In general, the magnetization is not aligned paral- 
lel to the z-axis. We therefore apply a coordinate 
transformation that rotates our system to align its 
z-axis parallel to the magnetization. The calcula- 
tions are much easier and also more convenient in 
the new system referred to as E'. 

2. After this we write down the equation of motion 
of the single magnon Green function. To solve this 
equation it is necessary to 

3. decouple higher operator combinations. This ap- 
pears to be straightforward for the exchange term 

as long as one works in the rotated system E'. 
We will perform the usual RPA decoupling! here. 

4. The situation is more complex for the anisotropy 
term J2J). Here we will develop a new decou- 
pling scheme following the ideas of the Callen 



decoupling^. However the original Callen decou- 
pling is not applicable, since the total spin is not 
a conserved quantity in our model. Therefore the 
decoupling has to be generalized. 

5. There is a special rotation angle 9 in our model: 
If the coordinate system is rotated by this angle 
E — > E and the decoupling procedure is applied, 
the total magnetization Sa commutes with the 

Hamiltonian It is easy to show that 9 is there- 
fore the direction of the magnetization. Now the 
condition [J^. Sa, H]- — gives an explicit ex- 
pression for the magnetization angle. 

6. Using the decouplings as well as the commutation 
property in the primed system we can solve the 
equation of motion and finally obtain the single 
magnon Green function as well as the norm of the 
magnetization (S z >). Therewith the problem will 
be solved. 

7. One can further show, that the effect of the 
anisotropy can be interpreted instructively as an 
effective "anisotropy field". We will calculate the 
components of this field. 

Lets now follow this program in more detail. The rota- 
tion of the coordinate system is described by E' = ME, 
where M is a rotation matrix. Due to the symmetry we 
may confine the rotation to the zcc-plane without loss of 
generality. This means that y' — y and that the polar 
angle 9 fully characterizes the rotation. 

/ cos6> -sin.0 \ 
M = 1 (4) 
\ sin 9 cos 9 ) 

The z'-axis of the new system E' is set to be parallel to 
the magnetization direction. This gives: 

(S*0 = (S y ,) = 0. (5) 

The magnetizations in the fixed system E can now be 
read off from Eq. 10} : 

(S x )=sm9(S z ,), 

(S z ) = cos 9(S Z ,). (6) 

(S z ) is the magnetization component normal to the film 
plane while (S x ) denotes the component parallel to the 
film plane. (S Z '), consequently, is the total magnetiza- 
tion. Of course, the angle 9 is a priori unknown. 

Next we want to write down the equation of motion of 
the single magnon Green function G'^E) = {{Sf; Sj')) 
defined in the new system. Applying the transformation 
Q to the Hamiltonian one readily finds: 

EG'^E) = ([St',Sr']_) + ({[S+',H].Sr')) 
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with 



(U) 

= 2{S z ,)S ij - 2j£ PiijiE) - T' Kj (E)) 
I 

+K 2 {cos 2 9- isin 2 9) T%(E), 
+ (B x0 sin 9 + B zQ cos 9)G' l0 {E) 

- (B x0 cos 6 - B zQ sin 9) ( (Sf ; 

+2K 2 cos 9 sin 6{{Sf;Sr')) 

+2K 2 cos9sm0((Sf 2 ;Sr')) 

-2K 2 cos 9 sin 9 ({S+'Sf; Sr')) 

- {(S+'Sf + SfS+'iSr'}). 



(7) 



The sum over I runs over the nearest neighbors of site 
i. The combination of trigonometric functions is a 
consequence of the rotation . 



pa/ 

ij 



To proceed the operator products in T' rnno and 
have to be decoupled. For the nonlocal products in 
the former we choose a symmetric RPA decoupling 
AB -> (A)B + A(B). Using Eq. 10 which is valid in the 
primed system we find 



jz/ c+/ Rp A I q \ c+i 



(8) 



This is the same result as found in the original RPA 
approach 7 for the " parallel limit" . Let us emphasize 
that this is only the case for the primed system £'. Thus 
the third point of our program is done. 

The crucial fourth step introduces a new approxi- 
mation scheme for the single ion anisotropy J2J. It will 
be used to decouple the higher Green function r"' as 
well as to treat operator combinations appearing in the 
commutator Sa, The latter is important for 

the fifth step of our calculation. 

1963 Callen introduced the Callen decoupling 9 which was 
intended as an improvement to the RPA decoupling 7 . 
The decoupling was performed at the Heisenberg ex- 
change term Later Anderson and Calient used 
this proposal to treat the local anisotropy term @ 
in the "parallel limit". Both approaches are based on 
the fact that the total magnetization commutes with 
the Hamiltonian, a condition that is not valid in our 
case. Therefore the procedure has to be rederived and 
generalized. We will not quite adopt the procedure as 
presented in Ref. 9 but rather use the main ideas. 
The approximation makes use of the operator identity 



The spin operators work on the same site i, the site index 
is dropped here for convenience. Now the "zero" is added 
to the components of the spin operator, S x , S y , S z : 

S(x,y,z) S(x,y,z) ~t~ 

a {x , y , z) (S(S + 1)-S 2 X - S 2 - Si) (10) 

It is important to note that theses relations are iden- 
tities for any (a^ j, z )) only for exact calculations. On 
the contrary the result of some standard approximation 
procedure (e.g. of a symmetric mean field decoupling) 
changes if one uses the right hand side of Eq. (fltjfl in- 
stead of the left hand side. The results do depend now 
on the prefactors {a( x ,y, z ))- It was the idea of Callen to 
use this degree of freedom to improve approximations. 
We will follow the proposal of Callen here and adjust 
the parameters in a way that interpolates between zero 
temperature and Curie temperature requirements. The 
explicit calculation is given in Appendix A. It gives: 

_ {S( Xy y,z)) 

a (x,y,z) — 



(11) 



Now we want to decouple operator combinations like 
S X S Z + S Z S X , S y S z + S z S y , or S x S y + S y S x which ap- 
pear in the equation of motion (0 in the higher Green 
function L°'. Thereto we replace the single operators 
by the right hand side of Eq. (|10f> and perform a sym- 
metric decoupling procedure at the resulting expressions. 
Since the prefactors a^ x y ^ are small quantities we ne- 
glect terms of the order a 2 . For example the result for 
the operator combination S y S z + S z S y is thus given by: 



SyS z + S Z Sy 



A.C. 



2{S V )S Z + 2{S z )S y (12) 
ASy) 



,{Sz)_ 

'2S 2 



(^{S X S Z ~\~ S z S x ) S x 

'(SySz + S Z S y )Sy + 2(S 2 )S Z ^J 

(^(S x S y + SyS X }S X 



+2{S 2 )S v + (S z S y + S y S z )S l 

For the other operator combinations analog expressions 
are found. 

By virtue of relation JSJ the result can be simplified if the 
decoupling is performed in the rotated coordinate system 
We show in Appendix B that the following relations 
are fulfilled in £': 

(S a >S V +S V S a >) = 0, (13) 

where a' and b' are two different subscripts out of 
(x',y',z ! ). This, together with Eq. JSJ, finally gives the 
decoupling within the parallel system E': 

S X 'S Z ' + S Z 'S X > — -> 



Sy' S z 



S Z 'Sy' 



., . , i , (- C y) 
's 2 ' 



2(S Z .) 1 - 



• Sy' 



S(S + 1) — s 2 — s 2 — S 2 Z 







(9) 



S X 'S Z , + S Z 'S X , ^ 2{S Z ') (l - ■ S xl (14) 
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Now the operator combination that appears in the Green 
function r°' in the equation of motion J7J can be decou- 
pled. Using Eq. 114(1 as well as the identity @ one finds: 



S+iS z ' + S z >S+> 



A.C. 



2(S Z ,)C[S. 



with 



C[ = 1-^(S(S + 1)-(S 2 Z ,)). 



(15) 



(16) 



Using the decoupling procedures discussed up to 
now the higher Green functions T' mno and Ffj in the 
equation of motion J7J) can be treated. To treat the 
other four terms we have to address the fifth point of 
our program. Hence we will show in the following that 
for a certain angle 9 the total magnetization J^. S? 
commutes with the Hamiltonian . Furthermore an 
explicit expression for the magnetization angle will be 
obtained. 

Applying the rotation (@J to the Hamiltonian © and 
using the abbreviation 71 = sin 9 and 72 = cos 9 we find: 



i 



Now the last operator product is decoupled according to 
Eq. ((HI)- This gives: 



S4 , H 



( 71-Bzo - 72-Bzo 



(17) 



+2K 2lll2 {S z ,) (1- 



S 2 



■iSV 



Thus, in the framework of our approximation, the total 
magnetization indeed commutes with the Hamiltonian if 
the term in brackets on the right hand side is zero. This 
has important consequences: All expectation values and 
Green functions in the rotated coordinate system £ that 
do not conserve the spin are zero in the framework of 
our theory. This can be seen using, e.g., the Lehmann 
representation of the Green function. In particular this 
applies for (S +/ ) and (S 1 -,). Therefore Eq. © holds in 
this coordinate system which is thus found to be equiv- 
alent to the system £': £ = Hence from Eq. ijT%|) 
follows simple condition for the magnetization angle 9: 



= sin 9B z0 — cos 9B xQ 

+2K 2 sin9cos9{S z ,)C[ 



(18) 



The expectation value (S 2 ) is already evaluated here 
using the property of spin conservation in the primed 
system. Eq. (|18H is our first important result. 



Having calculated the magnetization angle now the 
norm of the magnetization {S z i ) has to be derived. This 
turns out to be a straightforward task. Due to the 
property of spin conservation in the primed system the 
Green functions in the last four lines of the equation of 
motion Q are identical to zero. The remaining higher 
Green functions T' mno and can be decoupled by Eqs. 
(JHJ) and (|14fl . Thus the equation of motion can be 
solved after Fourier transformation. We finally obtain 
the single magnon Green function G'^E), 



with 



2(S Z ,) 
E-E' 



= 2(&/>J(p-7q)+S 



B = 



B x q sin 9 + B z0 cos 9 
1 , 1 



sm 2 9)2(S z ,)C 1 . (19) 



The term p denotes the coordination number, while 7 q is 
a structural factor due to the Fourier transformation of 
the Heisenberg exchange term. For the quadratic lattice 
chosen here it is given by: 

7 q = 2 (cos aq x + cos aq y ) , 

where a is the lattice constant. The trigonometric 
functions in Eq. I(19|) are obviously a consequence of the 
rotation. Knowing the Green function G'^E) one can 
calculate the desired expectation values in the primed 
system, i.e. the total magnetization (S Z ') and (S^,) 
by a standard text book procedure^ finally ending 
up with a self consistent system of equations. Before 
discussing the results of our theory in more detail 
we want to offer an instructive interpretation of the 
work of the anisotropy in the framework of our approx- 
imation. This will be the last point of our theory section. 

The abbreviation B in Eq. (|19l) has the same ef- 
fect on the Green function as an external field aligned 
parallel to the magnetization. Combining the expression 
for B with the magnetization angle, we may write down 
the components of the effective field: 

B x = Bsin9 

= B x0 - K 2 {S X ) sin 2 9C[ 

= B x q + B xa 

B 7 = Bcos9 



= B z0 + 2K 2 (S z )(l - isin 2 ff)C{ 
= B z0 + B za . 



(20) 



Obviously, the effective field may be written as a sum 
of the external field and an "anisotropy field" B a = 
(B xa , 0, B za ). The anisotropy acts exactly like this field 
as far as the magnetization and the magnon energies i? q 
are concerned. 
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FIG. 1: The x- and the z-component of the magnetization 
(S x ), {S z ) as a function of the external field calculated with 
our RPA+A.C. approach (solid line), the approximation pro- 
posed in Ref. 8 (dashed line) and with a mean field de- 
coupling of the anisotropy term (dotted line) in comparison 
with the QMC results from Ref. 6 (symbols). Parameters: 
K 2 = +0.02J and k B T = 2J « 0.32fc s T c , 5 = 2. 



In the next section we will present the results of our the- 
ory and compare them with QMC data and other ap- 
proximations. 



III. RESULTS AND DISCUSSION 

We start our discussion with a comparison of our re- 
sults to QMC data of Ref. 6, which are free of systematic 
errors. Fig. ^ shows the results for the field induced re- 
orientation transition at finite temperatures. A positive 
anisotropy constant (K^ > 0) and an external field paral- 
lel to the film plane (9b — 90°) are considered. The ex- 
ternal field is applied perpendicular to the easy direction 
of the magnetization, a situation representing a severe 
test for our theory. The components of the magnetiza- 
tion (5a;) and (S z ) are shown as functions of the external 
field. For zero field the magnetization is aligned parallel 
to the easy axis. It is fully rotated into the film plane at 
the reorientation field B t q. We display the results of our 
calculations (RPA+A.C. - solid lines) and the QMC re- 
sults of Ref. 6 (symbols). Additionally the results of two 
other theories are shown for comparison: For the dotted 
line the anisotropy term (0) is treated by simple mean 
field decoupling: 

5?5f ^ 2(S z )St (21) 

The dashed line shows the proposal of Ref. 8. Here the 
operator combination S^ Sf+S*S^ is decoupled as in the 
parallel limit treated in the original paper of Anderson 




FIG. 2: The z-component (5 Z ) (solid lines) and the x- 
component (S x ) (dashed line) as a function of the external 
field B x o applied within the film plane. Further parameters: 
K 2 = 0.01J, 5 = 2, k B T c equals 



and Calient: 

5+5* + 5?5+ - 2(5.) - ^(S(S + 1) - <5 2 2 ))) . 

The exchange term in all model calculations is treated 
by an RPA decoupling. 

The results of our theory (solid lines) are in excellent 
agreement with the QMC data. We achieved even quan- 
titative agreement for all magnetization angles 9. The 
quality of the approximation indeed turns out to be the 
same for all angles 9, which was the aim of this paper. 
Our approach is clearly superior to the approximation 
proposed in Ref. 8, to the mean field decoupling, and to 
all other approximations shown in Fig. 11 of Ref. 6. Fig. 
n visualizes our main result, namely that we succeeded 
to develop a theory for the extended Heisenberg model 
© that is as accurate as the RPA for the model Q. In 
the following we will discuss some additional features of 
the reorientation transition. 

In Fig. |2 the temperature dependence of the transition is 
analyzed. Again, the components of the magnetization 
are displayed as a function of the external field, which 
is applied in the film plane, perpendicular to the easy 
direction. The calculations were performed for three dif- 
ferent temperatures. Since the system is not saturated 
at finite temperatures, the total magnetization increases 
with the external field. This is seen best after the re- 
orientation (Bq > Ba r ), where only one component of 
the magnetization ({S x }) is present. For higher temper- 
atures the transition as a function of the external field 
becomes sharper. The reorientation field Bq t decreases 
faster with temperature than the zero field magnetiza- 
tion, reflecting the fact that the anisotropy becomes less 
important at higher temperatures. Another interesting 
feature is that the x— component increases linearly with 
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FIG. 3: The magnetization, the reorientation field Bq t (dots) 
and the anisotropy field B a for 9 = 0° (dashed line) and 
6 = 90° (dotted line) as a function of temperature. All quan- 
tities are scaled to their zero temperature value. The other 
parameters are as in Fig. |21 For Comparison scaled magneti- 
zation curves are added (solid line). 



FIG. 4: The z— and a;— component of the magnetization as a 
function of the external field. Results for positive K2 = 0.01J 
(solid lines) and negative Ki = — 0.01J (dashed lines) are 
shown, T = 0.7T C . The small arrows highlight the position 
of the reorientation fields B^ r for positive and B^ r for nega- 
tive anisotropy. The inset shows the magnetization curves for 
Bo = 0. The other parameters are as in Fig. |5] 



the external field until the reorientation field is reached. 
This holds for all temperatures and is qualitatively dif- 
ferent from the approximation proposed in Rcf. 8. Qual- 
itatively, this feature is also found in mean-field theory 
as can be seen in Fig. 1 of Ref. 14. 

The observed behavior follows directly from Eq. (|18fl . 
Since the external field is applied parallel to the film plane 
one finds for the x-component of the magnetization (S x ): 



sin8{S z ,) 

B x o 
2K 2 C[(T)- 



(23) 



Since C[ (Ea ll9|) increases with temperature, the slope 
of (5a;) is steeper for higher temperatures. Additionally, 
Eq. I|23|) determines the reorientation field Bo r . We find: 



B 0r (T) = 2K 2 (S Z ')(T) C[(T). 



(24) 



The fast decay of the reorientation field with tempera- 
ture as compared to the magnetization is also due to the 
temperature dependence of C' x . 

This can also be seen in Fig. |3 where we considered 
the temperature dependence of the system in detail. We 
plotted the norm of the components of the anisotropy 
field l|2l)[) as well the reorientation field (circles) and the 
magnetization (solid line) as a function of temperature. 
All quantities are scaled to their zero temperature value. 
The anisotropy fields are plotted at their maxima, i.e. at 
9 = 0° for B az (dashed line) and at 6 = 90° for B ax (dot- 
ted line) . The temperature dependence of the anisotropy 
fields (|20|l as well as of the reorientation field l|24|) are de- 
termined by the factor (S Z >)(T) C[{T). Thus this quan- 
tities have nearly the same temperature dependence and 
their slopes are steeper than that of the magnetization 



{S Z >}(T) alone. 

Very similar results are found for the easy plane case 
(K2 < 0, Bq||z). In Fig. 0] we compare both cases of a 
reorientation transition. Solid lines show the transition 
for an easy axis system, dashed lines denote the easy 
plane case. In the inset, the respective magnetization 
curves M(T) — (S Z >)(T) are plotted for zero external 
field. The Curie temperature and the magnetization at 
finite temperatures are somewhat smaller for the easy 
plane system. A reduced magnetization leads to a re- 
duced reorientation field (see Eq. (HU). This explains 
the differences between both cases concerning the reori- 
entation transition as seen in the main panel. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper we addressed the magnetic reorientation 
transition in a Heisenberg monolayer as a function of 
the external field and temperature. The basis of our 
approach is a transformation of the Hamiltonian into a 
coordinate system £' (with the z 1 axis parallel to the 
magnetization) as well as a generalized Anderson- Callen 
decoupling procedure. Compared to the bare Heisenberg 
Hamiltonian Q), the problem is more complicated, 
since the total spin is not conserved. However, this 
complication turns out to be less serious, as it can be 
shown that the total spin is a conserved quantity in 
the framework of our approximation, if a appropriate 
quantization axis is chosen. This fact can be used to 
calculate the magnetization angle as well as to solve 
the equation of motion for the single magnon Green 
function. It was further shown that the effect of the 
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anisotropy can be described by an effective "anisotropy 
field". 

Our results show a strikingly quantitative agreement 
with the QMC data of Ref. 6 yielding a significant im- 
provement over all other decoupling schemes discussed 
so far (see e.g. Rcf. 6). The main practical virtue of the 
new approach is that calculations can be performed as 
accurate as with QMC but much faster. 
The theory can be generalized to a multilayer system 
and can thus be used for cases where QMC calculations 
are not feasible any more (e.g. thicker films). It should 
therefore be used to analyze the magnetic reorientation 
transition as a function of the film thickness as found in 
many transition metal films (see e.g. Ref. 4). 
Due to its accuracy and convenience the theory shall 
further be used for a quantitative analysis of ferromag- 
netic resonance (FMR) experiments^i^. The decisive 
feature for the interpretation of a FMR experiment is 
the dependence of the q — spin wave mode -E q= o on 
the external field Bo- The function i? q= o(Bo) can be 
easily calculated in our theory for any direction of the 
external field. This opens the possibility to extract the 
microscopic anisotropy constant K2 directly from FMR 
experiments. 
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If a z is set to zero, the left hand side and the right hand 
side of Eq. (|25l) are approximated on the same level, 
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placed by their zero temperature value. However one can 
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On the other hand for high temperatures T > T c the left 
hand side of Eq. I|25(l has to vanish. This can be assured 
by the choice 

a,(T>T c )~coB0<S , J ,/}. (30) 

Combining the settings (|28|) and i|3U|l one ends up with 
Eq. (JTTJ: 

<_ co S 9(S z ,) _ (S z ) 
z ~ 2S 2 " 2S 2 ■ [ ' 

Analog calculations lead to the prefactors a x and a y . 
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Appendix A 

In this Appendix we derive the prefactors of Eq. (|10fl . 
We follow the philosophy of Callens papes^ and calculate 
the prefactors as an interpolation between low and high 
temperatures. Lets start with the former limit T ~ 0: 
Starting point is Eq. (|10fl . We will consider expectation 
values instead of operators and transform the resulting 
expression as: 

{S z ) = {S z )+a z {S(S + l)-S 2 x -S 2 y -S 2 ) 
= (S z ) + a z (S(S + 1) - Si, - S 2 y , - Si) 

= (s z ) 

+a z (S(S + 1) - S z , - S-,S +I - Si) (25) 

The primed terms are quantities of the rotated system 
which is aligned parallel to the magnetization. Now the 
expectation values of the right hand side are approxi- 
mated by their zero temperature values 

(S z -) S 

(Si) ^ S 2 (26) 

This gives 

(S z )ttScos6 + a z (-S-,S +l ) (27) 



Appendix B 

Here we want to derive the relation (|13|) . Starting point 
is Eq. JSJ). First we want to calculate the expectation 
value (S x 'Sy, + S y >S x >). Using the decoupling l(T2l to- 
gether with Eq. JSJ) one finds 

(S x ,S y , + S y ,S x ,) -> (S x ,)A+{S y ,)B 

= (32) 

The terms A and B are given by the decoupling proce- 
dure l|12l) . This is one of three equations that have to 
be derived to prove relation l|13fl. Next we want to treat 
(S X ,S Z , + S z ,S x i). Using the decoupling rule l|12|) as well 
as the result (|32H one finds: 

S X 'S Z , + S Z 'S X , > 2(S z i)S x , 

+ \(s x ,s z , + s z ,s x ,)s z ,) 

Thus it follows for the expectation value 

(S X ,S Z , +S Z ,S X ,) -» 2(S Z ,)(S X ,) 

-2^l(2(Sl,)(S x ,) 

+ l(s x ,s z , + s z ,s x ,) {Sz ,)) 



therefore: 


and 




2S 2 



(S X ,S Z , + S Z ,S X ,) 1 + 



— (S X 'S Z ' + S z 'S x i) 



((S X ,S Z , + S Z ,S X ,)(S Z ,)) 



2S 2 



(33) 



The equation 



(S y ,S z , +S z ,S y ,} =0 



(34) 



is derived in a analog way. Eqs. (J2HJ), and i|54"|) 

prove relation (|13l) . 
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